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We develop a general formalism for the parameter-space metric of the multi-detector J^-statistic, 
which is a matched-filtering detection statistic for continuous gravitational waves. We find that 
there exists a whole family of jT-statistic metrics, parametrized by the (unknown) amplitude pa- 
rameters of the gravitational wave. The multi-detector metric is shown to be expressible in terms of 
noise- weighted averages of single-detector contributions, which implies that the number of templates 
required to cover the parameter space does not scale with the number of detectors. Contrary to 
using a longer observation time, combining detectors of similar sensitivity is therefore the computa- 
tionally cheapest way to improve the sensitivity of coherent wide-parameter searches for continuous 
gravitational waves. 

We explicitly compute the .7^-statistic metric family for signals from isolated spinning neutron 
stars, and we numerically evaluate the quality of different metric approximations in a Monte-Carlo 
study. The metric predictions are tested against the measured mismatches and we identify regimes 
in which the local metric is no longer a good description of the parameter-space structure. 



I. INTRODUCTION 

Continuous gravitational waves (GWs), which would 
be emitted, for example, by spinning non-axisymmetric 
neutron stars, or by solar-mass binary systems, are gen- 
erally expected to be so weak that they will be buried 
several orders of magnitude below the noise of even the 
most sensitive detectors. The detection of such signals 
therefore requires the exact knowledge of their waveform, 
in order to be able to coherently correlate the data with 
the expected signal by matched filtering. 

In a wide-parameter search for unknown sources, how- 
ever, we typically only know the family of possible wave- 
forms (or an approximation thereof), parametrized by 
unknown signal parameters (such as the frequency or sky 
position of the source). The corresponding parameter 
space needs to be covered by a finite number of "tem- 
plates" for which a search will be performed. These tem- 
plates must be placed densely enough, so that for any 
possible signal, no more than a certain fraction of the 
signal-to-noise ratio (SNR) is lost at the closest tem- 
plate. On the other hand, coherently correlating the 
data with every templates is computationally expensive 
and increases the number of statistical false-alarm candi- 
dates. Therefore an optimal covering is desirable, which 
minimizes the number of templates but guarantees the 
required "minimal match" . 

In order to solve this covering problem, it is essential 
to understand the underlying parameter-space structure. 
Most studies on the construction of optimal template 
banks have been performed in the context of binary- 
inspiral searches. It was realized early on that a geo- 
metric approach to this problem is the most natural, in 
particular the introduction of a metric on the parame- 
ter space by Balasubramanian et al. and Owen 0, 
building on the earlier concept of the "fitting factor" in- 
troduced by Apostolatos Note, however, that this 
definition of the metric differs subtly from the "canoni- 
cal" definition used in the present work (and also in 3), 
which is derived directly from the detection statistic (see 



appendix IXI for more details). The canonical definition of 
the metric assigns the relative loss of SNR due to an offset 
in signal parameters as an invariant "distance" measure, 
which can locally be expressed as a metric tensor. The 
metric is closely related to the well-known concept of the 
"Fisher information matrix" , which quantifies the statis- 
tical errors in the parameter estimation of signals: the 
(canonical) metric is identical to the normalized Fisher 
matrix, even though it describes conceptually rather dif- 
ferent aspects of the detection statistic. 

A somewhat related question is the global parameter- 
space structure, which was studied in Q for the case of 
isolated neutron-star signals. This study found that the 
global structure (the "circles in the sky" ) deviates signifi- 
cantly from the local metric picture. The global structure 
is relevant, for example, for deciding whether different 
detection candidates are consistent with the same signal, 
i.e. whether they are "coincident candidates". Obvi- 
ously, the metric description is the local approximation 
to this global parameter-space structure. 

In this paper we consider gravitational-wave signals 
that are nearly monochromatic and sinusoidal in the 
frame of the GW source, and which are of long dura- 
tion (i.e. typically longer than the observation time T). 
This class of signals is usually referred to as "continu- 
ous waves" , and the prime examples are GWs from non- 
axisymmetric spinning neutron stars (e.g. see for a 
review) and stellar-mass binary systems in the LISA fre- 
quency band (e.g. see 3). 

The phase of the signal received at the detector is 
Doppler modulated by the rotation and orbital motion 
of the detector. The observed phase therefore depends 
not only on the intrinsic frequency evolution of the sig- 
nal, but also on its sky position. In addition to the phase 
modulation, there is a time-varying amplitude modula- 
tion, due to the rotating antenna pattern of the detec- 
tor. This amplitude modulation depends on the polar- 
ization angle tp and the polarization amplitudes ^+ and 
Aye of the GW. However, as shown by Jaranowski et al. 
[3|, these unknown parameters (together with the initial 
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phase (f>o), can be eUminated by analytically maximizing 
the detection statistic. The resulting reduced parameter 
space includes only the parameters affecting the time evo- 
lution of the signal phase, which we will refer to as the 
"Doppler parameters" . This amplitude-maximized de- 
tection statistic is generally known as the "JF-statistic" , 
which has been used in several searches for continuous 
GWs from spinning neutron stars (e.g. (8 , J3_, JJ)] ) . After 
two earlier (partly successful) attempts to generalize the 
jF-statistic to a coherent network of detectors 0, Q , this 
problem was fully solved more recently by Cutler and 
Schutz 

Somewhat surprisingly, however, there has not been 
much work on the metric of the .7^-statistic, neither in the 
single- nor the multi-detector case: the single-detector T- 
statistic metric was derived on a formal level by Krolak 
et al. 0, but was not evaluated explicitly or studied 
further. A single-detector JF-statistic metric was used 
(without giving any details) in to numerically esti- 
mate the number of templates in galactic-bina ry s earches 
with LISA. An earlier study by Brady et al. [13 of the 
metric for isolated neutron-star signals introduced a met- 
ric approximation based only on the phase modulation 
of the signal and neglecting the amplitude modulation. 
We will refer to this approximation as the "phase met- 
ric". This metric has a simpler structure than the full 
^-metric, and it can be computed analytically if one 
assumes a circular orbital motion. This is the only type 
of continuous GW metric that is currently implemented 
in LAL/LALApps JJ]- 

As we will see in this study, the amplitude modulation 
cannot always be neglected, but "on average" the phase 
metric seems to be a good approximation, and its qual- 
ity improves with longer observation times and with the 
number of detectors. With the recent multi-detector gen- 
eralization of the J^-statistic formalism 11] and its subse- 
quent implementation into LAL/LALApps, the need to 
understand the multi- detector JF-statistic metric has be- 
come more urgent. The most important question in this 
context is whether the metric resolution increases with 
the number of detectors, i.e. whether a denser covering 
of the parameter space is required, which would increase 
the computational cost. 

The main result of this work is to show that the metric 
resolution does not scale with the number of detectors. 
Therefore, sensitivity can be gained at the cost of only a 
linear increase in the required computing power (as the 
signal has to be correlated with the data stream from 
each detector). This has to be contrasted with the (at 
least) 0(T^) scaling l|117|l of the number of templates 
with observation time T, in the case isolated neutron- 
star signals with one spindown. 

In order to improve the sensitivity of a coherent search 
for continuous gravitational waves, increasing the num- 
ber of similar-sensitivity detectors is therefore computa- 
tionally much cheaper than to increase the observation 
time. 

The plan of this paper is as follows: in Sect. ^ we in- 



troduce the formalism and notation of the multi-detector 
J^-statistic, following and H. In Sect. UTTl we derive 
the JF-statistic metric family for high-frequency signals 
(relevant for ground-based detector networks). We com- 
pute the extremal range of this metric family, its aver- 
age metric, and the long-duration limit, in which the J-- 
metric family reduces to a simple "orbital metric" . In 
Sect. lIVl we apply this framework to the special case of 
GWs from isolated spinning neutron stars, and we eval- 
uate the quality of the metric predictions (and different 
approximations) by comparing them against measured 
mismatches in a Monte-Carlo study. The main results 
are summarized in Sect.0 Appendix 1X1 presents an al- 
ternative, more elegant derivation of the J^-statistic met- 
ric, and Appendix IbI gives the general expressions for the 
J^-metric, which would be valid also for low-frequency 
signals relevant for LISA. 



II. THE MULTI-DETECTOR JP^-STATISTIC 

A. General definitions 

In this section we introduce the formalism and notation 
of the JF-statistic, a matched-filtering detection statistic 
for continuous gravitational waves, which was first intro- 
duced by Jaranowski et al. , and subsequently general- 
ized to the multi-detector case by Cutler and Schutz . 
As shown in 7], the dimensionless strain signal s'^{t) of a 
continuous gravitational wave at detector X can be rep- 
resented in the form 

4 

s^(0 = ^^^/i^(t), (1) 

in terms of four signal-amplitudes ,4'', which are inde- 
pendent of the detector X, and the detector-dependent 
basis waveforms h^{t). The four amplitudes A'^ can be 
expressed in terms of two polarization amplitudes A^, 
Ax, the initial phase 00 in the solar-system barycenter 
(SSB) at a reference time Trof, and the polarization an- 
gle ip of the wave frame with respect to the equatorial 
coordinate system, namely 

= cos 00 cos2tl; — Ax sin^o sm2ip , 

A? = cos 00 sin 2?/; + Ax sin 00 cos 2 , 

A? = —Aj^ sin 00 cos 2-0 — Ax cos0o sin 2?/;, 

A^ = — sin 00 sin2'0 + y4x cos0o cos2?/'. 

We can further relate the two polarization amplitudes A_|_ 
and Ax to the overall amplitude ho and the inclination 
angle l of the quadrupole rotation axis with respect to 
the line of sight, namely 

A+ — \:ho (1 -|- cos^ i) , Ax = ho cos i . (3) 
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The four basis waveforms h^{t) can be written as 

hf{t) ^a^{t) cos ^^{t), hf{t) ^b^{t) cos (l)^{t), 
hf (t) = (t) sin 0^ (t) , hf{t)^b^ (t) sin 0^ (t) , 

where a^{t) and b^{t) are the antenna-pattern func- 
tions (see Eqs.(12,13) of 0), and (jP^{t) is the signal 
phase at the detector X. The antenna-pattern functions 
a?^{t), b^{t) depend on the sky position n of the GW 
source, and on the location and orientation of the de- 
tector X. The phase (ly^{t) also depends on the intrinsic 
phase parameters, cu say, of the signal. In the case of con- 
tinuous waves from isolated neutron stars, lo would only 
consist of the s + 1 spin parameters, i.e. ui = {/'''■'Ifc^o' 
where /''^-' is the k-th time-derivative of the intrinsic sig- 
nal frequency in the SSB (see Sect. IIV|l. In the case of 
signals from spinning neutron stars in a binary system, 
uj would also contain the binary orbital parameters. We 
can summarize these dependencies as 

a^ = a^{t;n), b^^b^{t;n), (jp^ ^ (j)^{t;n,Lj) . (5) 

In the following we denote the set of "Doppler parame- 
ters" (i.e. the parameters affecting the time evolution of 
the phase) by A = {n,Lu}, as opposed to the four "ampli- 
tude parameters" {.4}^ — . Note that in the litera- 
ture the Doppler parameters are sometimes referred to as 
"intrinsic" , and the amplitude parameters as "extrinsic" , 
but we will not use this convention here. 

Using the multi-detector notation of 0, ^ '^^ write 
vectors in "detector-space" in boldface, i.e. {s}^ = s^, 
and so the signal model ^ can be written as 

sit;A,X) ^A^h^it;X), (6) 

with implicit summation over repeated amplitude in- 
dices, ^ G {1,2,3,4}. We assume the data x'^{t) of 
detector X contains a signal with parameters {A, A} in 
addition to Gaussian stationary noise n^{t), i.e. 

x{t) ^ n{t) + s{t; A, X) . (7) 

In general, the noise contributions n'^{t) of the different 
detectors could be correlated (which might be relevant 
for the two LIGO detectors in Hanford, or for LISA), 
and so we define the (double-sided) noise-density matrix 
as 

/OC 
^^^(r)e-^^-f^dr, (8) 
-oo 

in terms of the correlation functions 

(r) = £; [n^ (t + r) (t )] . (9) 

The corresponding multi-detector scalar product is de- 
fined (in analogy to O^) as 

/OC 
^''{f)s^Uf)y^*U)df, (10) 
-oo 



where x{f) denotes the Fourier-transformed of x{t). We 
use implicit summation over repeated detector indices, 
and the inverse noise matrix is defined by S^y S^"^ = 5\. 
In the case of uncorrelated noise, where = S-^^ , 
the scalar product simplifies to 

ix\y) = J2{xV), (11) 

X 

in terms of the usual single-detector scalar product 

With the definition (|10|l of the multi-detector scalar prod- 
uct, the likelihood function for Gaussian stationary noise 
n{t) can be written as 

P(n(i)|^^Y) =A:e-5("l"), (13) 

where A: is a normalization constant. Using Q we can 
express the likelihood of observing data x{t) in the pres- 
ence of a signal s{t; A, A) as 

P{x\A, A, S^^) = k e-5(^N) e(x|^)-i(«|«) (^4) 

Using a standard frequentist approach, the optimal de- 
tection statistic to decide between the hypothesis of a 
signal s — s{A, A) being present and no signal, s = 0, is 
given by the likelihood ratio A{x;s), defined as 

\nA{x;s) ^ In = {x\s) ^{s\s) . (15) 

Substituting the signal model , the log- likelihood ratio 
is found as 

\iiA{x;A,X) = A''x^-^M'M^.,A'', (16) 
where we defined 

x^(A) = ix\h^) , (17) 
M^^iX) EE {h^\K) . (18) 

We see that the likelihood ratio (|15|l can be maximized 
analytically with respect to the unknown amplitudes A'^ , 
and the resulting detection statistic for the Doppler pa- 
rameters A is the so-called -statistic, namely 

2T{x-X) = x^,M^''' x^, (19) 

where M^"' = {M-^Y" , i.e. Mf.o^M"'' = S^. Note 
that the four (multi-detector) wave- functions h^{t) form 
a basis of the signal space, and M is the associated 
metric, which allows us to raise and lower amplitude 
indices. In particular, we could define a "dual" basis, 
h^{t) = Mt"' Kit), satisfying {hf'\K) = Jff, and the 
statistic H19|) could then be written even more compactly 
as 2J^(A) = x'^ x^. 
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B. ^-statistic of perfectly matched signals 

Let us consider the case where the target Doppler pa- 
rameters A are perfectly matched to the signal Ag, so the 
measured data would be x{t) = n{t) + s(i; A, A). In this 
case, the projections (fT7|) are 

(yl, A) = ( A) + {A, A) , (20) 

where we defined = {n\h^) and = {s\h^). Assum- 
ing Gaussian stationary noise, one can show that 

E [rif,] = , and E [n^ n^] = M^,^ , (21) 

where E[] is the expectation value. We further find 

E[Xfj\^s^, and E[Xf_,x^] = Mf_,y + s^j_Sy , (22) 

so the four random variables a;^ have mean and co- 
variances M.fiu- Using this together with p9|l . we find 
the expectation value of the jF-statistic as 

E[2T]=A + p^{0), (23) 

in terms of the "optimal" signal-to- noise ratio p(0), given 

by 

p2(0) = M^"" s, = A^" M^. A" = {s\s) . (24) 

Using (|22|l . it is straightforward to show that the 
quadratic form (|19|l can be diagonalized as the sum of 
four squares of Gaussian random variables with unit co- 
variance matrix. This implies that 2J- is distributed 
according to a (non-central) x^-distribution with four 
degrees of freedom, as previously shown by Cutler and 
Schutz [Tlj . The corresponding non-centrality parame- 
ter of this x^-distribution is p^(0). 

C. jT-statistic of mismatched signals 

In the case of unknown signal parameters, there will 
generally be an "offset" AA between the Doppler param- 
eters As of the signal and the target parameters A, i.e. 

A = As + AA. (25) 
In this case, the projections H17() can be expressed as 

{A, As ; A) = ( A) -I- ^" Uo,^ ( As ; A) , (26) 
where the 4x4 matrix 7?.q^ is defined as 

7^„^(As;A) = (/i„(As)|VW) ■ (27) 

This matrix is generally not symmetric, but evidently 
satisfies the symmetry relation TZafj.{K', A) = 7?.^q(A; As). 
Using this together with l|19(l . we can write the mis- 
matched (SNR)2 as 

p2(^,As;A) = E[2T]~A 

= ^"7^„^(As;A)A^'^''(A)7^;3,(As;A)^'3.(28) 



Assuming the target parameters A to be "close" (in a 
suitable sense) to the signal, we can Taylor-expand these 
matrices around the signal- location As, and keeping only 
terms up to second order, we obtain 

+^9y-M'^''(As)A^AA^ (29) 
7^^,(As;A) = X^,(As)+7^^..(As)Ay 

+^^M-J(^s)AA*AA^ (30) 

where indices i,j refer to the Doppler parameters A* 
(with automatic summation), and where we defined 

TZ^^i = {hf,\dih^) , (31) 
7^^^,,y■ = {hf,\dijh„) , (32) 

with di = d/dX'- and dij = d"^ /dX'-dV . Substituting 
these expansions into H28() . we find to second order 

p2(AA) = [M^, + L^,,AA' - AA'AA^] A" , 

.(33) 

where the first- and second-order coefficients explicitly 
read as 

U ^ n^ + '^Ui+M-d^M-^ -M, (34) 

-g,, = i(7^„+T^„) + iM.a,,M-l.>[ 

+M ■ d,M-^ ■ ^Uj + Ui ■ ■ ^TZj 
+n^-djM-^-M, (35) 

using matrix notation, and writing for the transpose 
of the amplitude indices. By applying derivatives to the 
identity • A^^^ = I, we find 

M-diM-^-M ^ -d.M, (36) 
M-d^jM-^-M ^ -^^JM + 2^^M■M-^■^■jM,i37) 

and using (|18|l and the product rule, we further obtain 

d,M = n^ + ^n,, (38) 

d,.jM = +^n,j +2h,j , (39) 
where we defined 

h^i,ij = {dihf,\djh^) . (40) 

For simplicity of notation, we assume the Doppler indices 
in H35|) to be implicitly symmetrized, i.e. Qij = Qji, as any 
non-symmetric parts will not contribute to the quadratic 
form H33|) . Using the identities H36|) - (|39|l . we can reduce 
the expansion coefficients (|34|) . H35|) to 

L, = 0, (41) 
Qrj = h,,-^n,-M-^-nj. (42) 

The exact vanishing of the first-order coefficient shows 
that the perfectly-matched case, A = As, is a local ex- 
tremum of the J^-statistic, as expected. The zeroth-order 
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term in (|33(l corresponds to the perfectly-matched case 
(|24|l . and so we arrive at 

p\AX) = p\0)-A^gf,,.,,A'' AA^AA^+0(AA3), (43) 

where Gf^uij can be written more exphcitly as 

Q^..^J = {d,h^,\d,K)-{h^\d,h^,)M^P {hp\d,K) . (44) 

As discussed in 4] , this matrix is directly related to the 
projected Fisher matrix F, namely 



(45) 



describing the statistical information on the Doppler pa- 
rameters A*, given the amplitude parameters A'^. 



III. THE jr.STATISTIC METRIC FAMILY 

A. General definitions 

The relative loss in expected JF-statistic due to an offset 
AA with respect to the signal location As defines a natural 
dimensionless "mismatch"^ namely 



mjr{A, As; AA) 



p^(0)-p^(AA) 
p2(0) 



(46) 



Using (|43() for the mismatched SNR, we can cast the !F- 
mismatcli in the form 



= gf^ {A, A) AA'AA^ + 0(AA^) , (47) 
where we defined the JF-statistic metric gf, as 



gijiA, A) 



A-g^jiX)-A ^ 
A-M{X)-A p2(o) 



(48) 



This expression is identical to that found previously 
(for the single-detector case) in 0, 0|, referred to as 
the "normalized projected Fisher matrix" . A more ele- 
gant method of obtaining this result by projecting the 
full parameter-space metric into the Doppler subspace is 
shown in appendix^ It is obvious from (|46|l that the 
overall signal amplitude ho cancels out, and that the mis- 
match mjr is therefore, contrary to the Fisher matrix F^ , 
independent of Hq. 

The Fisher matrix Tij characterizes the statistical un- 
certainty (due to the presence of noise) of the maximum- 
likelihood estimators for the signal parameters. In par- 
ticular, the inverse Fisher matrix gives lower bounds on 
the variances of the parameter estimators (the so-called 
Cramer-Rao bound), describing the typical fluctuations 
by which the location of the detection-statistic peak will 



^ See appendix ^ for a discussion of a slightly different mismatch 
definition sometimes found in the literature. 



vary in repeated experiments. The Fisher matrix can 
therefore be regarded as a measure of the best possible 
accuracy of parameter estimation. The concept of the 
metric gf^, on the other hand, describes the relative (i.e. 
SNR-indepcndent) "extent" of detection-statistic peaks 
for any single realization of the noise. This gives a mea- 
sure of the intrinsic parameter-space resolution, which is 
not related to the true location of any putative signal. 
However, while conceptually rather different, Eq. H48() 
shows that there exists a deep connection between the 
two concepts. 

The J^-metric (|48|l is not a unique metric on the 
Doppler-parameter space A of the JF-statistic, as it de- 
pends on the (generally unknown) signal amplitudes A^ . 
Expression H48I) therefore describes a whole family of 
metrics, corresponding to different A — const, subspaces 
of the full parameter space, which is of very limited direct 
use for the covering problem of the Doppler-parameter 
space. 

However, we can explicitly compute the possible range 
of mismatches for any given A and AA. For this, consider 
the extrema of tojf as a function of A, i.e. 







dm 



2gA 



AGA 



dA AMA {AMA)'^ 



2MA, 



(49) 



where we wrote G = ^^ AA'AA^, which is a 4x4 matrix 
in amplitude space. Equation is equivalent to 



(M-^-g) A ^ mjr{X, AX) A , 



(50) 



which determines the extremal values fhjr of the JF- 
mismatch as the eigenvalues of Ai~^-g. According to 
the Rayleigh principle, the minimum and maximum of 
the mismatch will be given respectively by the smallest 
and largest eigenvalues. As suggested in a more 

practical mismatch measure can be constructed from the 
mean of the eigenvalues, and we can define an "average" 
J^-metric as 



\Tr[M- 



(51) 



where the trace Tr refers to the amplitude indices. This 
average .F-metric, contrary to gf^, is independent of the 
amplitudes A^ , and is therefore of more practical inter- 
est as a metric on the Doppler-parameter space. The 
corresponding average J^-mismatch mjp is simply 



m^(A, AA) = gfjiX) AVAA- 



(52) 



We note that in the analogous case of binary-inspiral 
signals, where one can equally maximize the detection 
statistic over some of the signal parameters (referred to 
as "extrinsic" parameters) , the metric of the reduced pa- 
rameter space depends again on both extrinsic and intrin- 
sic parameters [iM Il9j . A common choice of mismatch 
metric for the template placement in this binary-inspiral 
context is to use the most conservative case (based on the 
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concept of the "minimal match" namely the worst- 
case mismatch 

m^'''^(A, AA) = max^ m^(y^, A; AA) . (53) 

The is referred to as the "minimax" prescription in |18|. 
Contrary to the average metric H52I) , this extremal metric 
cannot be expressed as a quadratic form in the Doppler 
separations AA, and so the corresponding iso-mismatch 
surfaces are not described by hyper-ellipsoids. 

As we will see in the following section, for the type 
of narrow-band continuous- wave signals considered here, 
the 4x4 matrix A4^^ -Q has only two independent eigen- 
values, corresponding directly to the maximum and min- 
imum possible J^-mismatches for any given A and AA. 



B. Narrow-band signals, uncorrelated noise 

In the following we restrict our analysis to continuous 
gravitational waves with a well-defined, slowly varying 
(intrinsic) frequency /. This assumption applies, for ex- 
ample, for GWs emitted from spinning non-axisymmetric 
neutron stars, and from stellar-mass binary systems (rel- 
evant for LISA). We assume the observation time T to 
be much longer than the GW period l//s, such that the 
number of cycles N is large, i.e. 



TV = /, T > 1 



(54) 



The phase (f>{t) of the signal will be dominated by 
the zeroth-order term 2Trfst, while the intrinsic fre- 
quency variability and the Doppler modulations are much 
smaller corrections. Assuming two such narrow-band sig- 
nals x{t) and y(t), we can approximate the scalar product 
(Cni) as 



(a;|y)«T^xj(/,)(x^yY) 



(55) 



in terms of the time-averaging operator (.), defined as 



(g) 



1 

T 



g{t) dt . 



(56) 



Note that the noise matrix S'^^ {fs) can be considered 
as a metric in "detector space", allowing us to lower 
and raise detector indices, e.g. we could write xx{t) = 
S^yifs)^^ (t): and the scalar product ifK^ would then 
read as {x\y) =T {x-^ yx)- 

For simplicity, we restrict our analysis to the (more 
common) case of uncorrelated detector noises, i.e. we 
assume S^y = S^^Sxy, which allows us to introduce 
"noise weights" wx as 



wx 



s- 



-, with S-' = J2Sx\ (57) 



such that J2 = 1- Using this, the scalar product H55|l 
can now be expressed as 



{x\y)^TS-^{xy)s, 



(58) 



where we introduced the noise-weighted averaging oper- 
ator defined as 



(59) 



The scalar products involved in the expression for 
the J^-statistic consist of products of slowly vary- 
ing antenna-pattern functions g{t) £ {o{t), b{t)}, 
and highly oscillatory, quasi-periodic functions 
p(<) G {sin^ 0(i), cos^ (/)(i), sin0(t) cos0(t)}. In the 
following we approximate the oscillatory functions as 
exactly periodic with period t — 1//, i.e. </> « 2t: f t, 
and so we can write the time-average (gp) as 



1 ^1 rtr^+i 



dt. 



(60) 



where t„ = n r, and N — T/t is the number of oscillation 
cycles. We can Taylor-expand g{t) in each of the periods 
as g{t) = g{tn)+g{tn) it - 1„) \g{tn){t - 1^)^ + and 
using the fact that 



(t-i„)"p(t)dt=(t"p). 



(61) 



where (.)^ denotes the time-average over one cycle oip(t), 
we obtain the expansion 



T 



« {p)r{9) + ^[9iT)-gm- 
For the averages over one cycle r we find: 



(62) 



(sm 

(i™ sin^ (/))^ - (i™ cos2 0)^ - (i™ su 



cos'-0)^ = |, (sin0 cos(/))^ = , 



(63) 



Under the above assumptions, we can therefore write 

(gp) « (.9> (P> + a/N) , (64) 

and so we can neglect higher-order contributions and 
keep only the zeroth-order term (p) (g) in expressions 
of this type. Using this approximation, the 4x4 matrix 
M^i„, defined in lfTH|) . is explicitly found as 



where C is the 2x2 matrix 



C = 



A C 
C B I ' 



(65) 



(66) 



in terms of the three independent antenna-pattern coef- 
ficients 



A^{a\, B^{b\, C^{ab\ 



S ' 



(67) 
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and we further define D = AB — . Inserting tlie ex- 
plicit expressions ||2J) for the amplitudes and using 
()65|l , we can write the optimal SNR (|24|l explicitly as 



P'(0) = ]^hlTS~^[axA + a2B + 2a:iC] 
in terms of the amplitudes (writing r/ = cos t) : 
ai(77» = i(l + 77^)2 cos^ 2-0 + ?7^sin^ 2-0, 

03(77, V") = 7(1 - 77^)^ sin 20 cos 20 . 



(68) 



(69) 



We see that the optimal SNR does not depend on the 
initial phase 0Oj and it scales linearly with the overall 
amplitude ft-o, and with the square-root of the observation 
time T . The dependence on the number of detectors will 
be discussed in the next section. 



C. Dependence on the number of detectors M 

A question of central importance is how the SNR and 
the metric resolution depend on the number J\f of co- 
herently combined detectors. The dependence of the op- 
timal SNR is very easy to see: in the explicit expres- 
sion H68|l . the antenna-pattern coefficients A, B, C are the 
noise- weighted averages (|67|l over detector-specific quan- 
tities, and the only scaling with the number of detec- 
tors Af therefore comes from the total inverse noise floor 
= "^x^- If assume, for simplicity, that all A/" 
detectors have a similar noise floor Sq, i.e. S^^ 
then the optimal SNR scales as 



p(0)cx 



ho 



'So 



VrJf. 



(70) 



Doubling the number of detectors (of similar noise floor) 
therefore has the same effect on the SNR as doubling the 
observation time T. 

It is not difficult to see from (|58|l and definitions H18|) 
and H44() . that both the numerator and denominator in 
(|48|l have the same scaling with S^^ (and therefore Af), 
which cancels out. To show this more clearly, we write 
the explicit expression found from Eq. H48|l for the multi- 
detector JF-metric, namely 



9'. 



A^ 



{d.ih^^d.jK)g - {hadiK)g{hahi3)g ^{hpdjK) 



A" 



AP {hpha)s 



It is evident from this expression that the JF-metric only 
depends on noise-weighted averages of single-detector 
contributions, but does not scale with Af. Note, however, 
that the multi-detector metric is not a simple average of 
single-detector metrics. 

Increasing the number of detectors therefore does not 
increase the metric resolution in parameter space. This 
is in strong contrast to the effect of increasing the obser- 
vation time T, in which case the metric resolution, and 



therefore the number of templates, grows at a high power 
of T (e.g. cx 0(T^), see (|117ll '). We can therefore gain 
SNR cx VJ7 with Af similar-sensitivity detectors, at the 
cost of "only" a linear increase cx A/" in the required com- 
puting power. Using a coherent multi-detector search is 
therefore the computationally cheapest way to increase 
SNR in a coherent wide-parameter search for continuous 
gravitational waves. 
As 



discussed in Sect. IIII Al the metric resolution 
must not be confused with the accuracy of parameter- 
estimation. The latter increases with SNR (and therefore 
also with the number of detectors), as described by the 
Fisher information matrix, while the former does not. It 
might still seem somewhat surprising that the additional 
"information" coming from the time delays between de- 
tectors does not result in a higher metric sky resolution. 
This can be understood in terms of the diffraction limit, 
which can be used to estimate the order of magnitude of 
the sky resolution from the ratio of the wavelength c// 
to the "baseline" VT, where V is the orbital velocity. 
The expected metric sky resolution is therefore of the 
order Aflo ^ c/{TfV), which is exactly the metric scal- 
ing found in Sect. IIV Bl It is therefore evident that for 
detector distances of the order ~ 1, 000 km, the "integra- 
tion baseline" V T dominates the sky resolution starting 
from observation times as short as T > 100 s, and no 
additional sky resolution can be gained from the baseline 
spanned by different detectors. 



D. The .F- metric family for high-frequency signals 

In order to explicitly calculate the J^-metric family 
((7T|l . we need to consider the derivatives dih^ of the basis 
functions namely (omitting detector indices) 



dih2 



diO cos ( 
dib cos (, 
diQ sin(/ 
dib sin d 



— a di(j) sin < 

— bdi(f) sinc^ 
-t- adicf) cos( 
-I- bdiS cos d 



(71) 



The antenna-pattern functions a{t),b{t) do not depend 
on any of the Doppler parameters A except for the sky po- 
sition ft. From their explicit expressions (cf. Eqs.(12,13) 
of 0) for ground-based interferometers one sees that 



0(1), d.Eb^b^0{l), 



(72) 



while the corresponding phase derivatives will typically 
be of order 



\dnm\ - |9t0 



Ml 

c 



27r/ 



2TTfTVlc - 



(73) 



Ignoring the constant term, which will not contribute to 
the metric, the second term will be much larger than 
unity if the number of cycles N = fT satisfies N ^ 
c/V ~ 10^. This will always be be true for high-frequency 
signals relevant for ground-based detectors, for which we 
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can therefore neglect the antenna-pattern derivatives in 
(|71|l with respect to the phase derivatives, i.e. 



di(t) > dia , dib , 



(74) 



This assumption might not hold in the case of low- 
frequency signals that would be more relevant for LISA, 
or for very short observation times. The correspond- 
ing calculations are somewhat more tedious, but lead to 
equivalent results and are presented in appendix^ Us- 
ing (|74|l . (|58|l and keeping only the leading-order terms 
in H64(l . we can approximate (|4()|l as 



' pi pi 





V 



\ 



I 

V 1-3 / 







(75) 



with the three independent components 



Pl^ = {a^dicj,d,4>)s, 



P2 = 

a b did) di 



s ' 



(76) 



In the same way, can be approximated as 



where 




(77) 



R} 



Hi 



R}'^ = {abd^ 
Using this and Ht)5|) . we obtain 



'>5- 



(78) 



1 



\ 


V Qi;. J 



Qij Qij 



with the three independent components 

DQlj = A{abd,(l))s{abdj<j))g + B{a^di(j))g 
-2C{a^d^(f))g{abdj(t))g, 



(79) 
(80) 



DQl = A{b^d,(j))s{9djdy)s + B{abd,(f)s{abdj^). 



-2C{abd^4>)s{b'^dj<j))g, 



(81) 



A{abdi(t))s{b^dj(j))g 



B {abdi4i)g (a^ dj <j)) r 



-C [{b^di(j))g{a'^dj(j))g + {abd^(j))g{abdj(j))g\ ,(82) 



assuming implicit index symmetrization in Combin- 
ing this with H75|l and 144() . we find the 4x4 amplitude 
matrix in the form 
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/ w} 

in? m? 

Q nn} m? 

V Q m? m? 



(83) 



in terms of the three independent mismatch components 

wr{\; AA) = m[^- AA'AA^ , (84) 

where r €E {1, 2, 3}, and the corresponding matrices are 

ml^{X) ^ Pl^ - Ql^ . (85) 

In analogy to H68|l . we can write the projected Fisher 
matrix (|45|l as 

= i h^TS'^ [ai mlj + a2 mf^ + 2a3 mf^] , (86) 

in terms of the amplitudes ar — ar{cos l,^) defined in 
(|69|l . This allows us to express the J^- metric family (|48f) 
in the more explicit form 



"2 mfj + 2a3 mfj 



gf. ( cos L,ip;X) = — — — — 

^ aiA + a2B + 2a3C 



(87) 



As discussed in Sect. IIII AI the extrema rhjr of the mis- 
match family mjr(^A) are given by the eigenvalues of 
Ai^^ ■ Q, namely 

= det [M-'^ -g -mjrl] 

= detM-^ det[g -mjrM] , (88) 

and therefore 

= (m^ - fh3rA){m'^ - fhjrB) - {m^ - rhjrCf . (89) 

We see that there are maximally two independent eigen- 
values, namely 



- max mill 



m% 



where 



= {2D)-^ [B 



m 



m 

2 



+ Ani^ 



2Cm^ 



(90) 

(91) 
(92) 



Note that m? > is necessary for the positivity of the 
JF- mismatch mjF, and is equivalently to the matrix Qij 
being positive definite, as can be seen from H83I) . The ex- 
tremal solutions H90II determine the maximum and min- 
imum possible mismatch, respectively, as well as the av- 
erage .F-mismatch m^r, for given signal location A and 
offset AA. The corresponding average JF-metric is found 
from (|52Jl and (jHU as 



gfj{X) = {2D)-^ + Am% -2Cm 



y J 



(93) 



Note that, contrary to the average mismatch mjF, the 

- mill I max 



extremal values cannot be written as quadratic 

forms in the Doppler offsets AA. The range of possible 
J^-mismatches for given A and AA can be characterized 
by an intrinsic "uncertainty" Atojt, which we define as 



Am^(A; AA) = ^ (m 



max 



= y'm^, - m2 , (94) 
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and we further introduce the relative uncertainty ajr as 

1/2 



o-jfIA; AA) = 



1 



(95) 



which is bounded in ajr e [0, 1]. The maximal mismatch 
fii^^^ can never be larger than twice the average mis- 
match TOjF, and in most cases it will be smaller. The 
average ^-metric gfj might therefore be quite a reliable 
mismatch measure in practice, which will be confirmed in 
the Monte-Carlo studies in Sect. II VI for isolated neutron- 
star signals. 



E. Long-duration limit: the "orbital metric" 

Consider the limit of very long observation times com- 
pared to a day r^, i.e. T 3> t^. We can always write the 
signal phase 4P^{t) at the detector X as 



(96) 



where (f>oTbit) is the signal phase modulated by the orbital 
motion, while A(p^{t) accounts for the additional diurnal 
phase modulation due to the spin of the earth, which is 
of order 



A(/)^(i) - 2TTfAT^{t) 



(97) 



where the time delay |Ar'^(t)| < 0.02 s is periodic over a 
day. We restrict ourselves to signals with slowly varying 
frequency over the timescale of a day (e.g. excluding 
neutron stars in close binary systems), for which we can 
assume 4'oThit) ~ '^i^ft- With these assumptions and 
(|64|l . we can approximate the typical contributions to 
the J^-metric in the following way: 



X 

E 



(98) 



According to our assumptions, the average {di4>orh) is at 
least linear in the observation time T, while the second 
term will be roughly constant on timescales longer than 
a day. In the limit T 3> we therefore find 

{a^d,(b)s "^^'^ A (a,0orb) . (99) 
In the same way, we can approximate in this limit 



(100) 



where we defined 



lorb 



(a,0orb), and 0"^= (a0orb9j0orb). (101) 



Introducing the "orbital metric" g™^ as 



(102) 



and the corresponding orbital mismatch morb, 

mo,b = <?°f AA^ AA^ , (103) 

we find the limiting T 3> approximations 

TO^ « A morb , "rn^ ~ B morb , m-"^ ~ C JTiorb • (104) 

The J^-mismatch range (|90|) therefore reduces to a single 
eigenvalue, namely 



~ mill max 

rrijr w morb , 



(105) 



and the multi-detector JF-metric family reduces to the 
orbital metric (|102|l . i.e. 

g,^(AA)'^«'''<?,f (A), (106) 

which is independent of the unknown amplitudes and of 
the number and position of detectors. 

Note that a different, but related metric approxima- 
tion is the "phase metric" , which neglects the amplitude 
modulation a^{t), b^{t), but retains the detector-specific 
phase modulation AtjP^lt), namely 



4^(A) ^ (^0^9,0^) - (9.0^) (9,0^) 



(107) 



While this is intrinsically a single-detector metric, one 
could formally generalize it by analogy with the J-"- 
metric, and simply replace the time averages (.) by noise- 
weighted multi-detector averages In this manner we 
could define an ad-hoc multi-detector phase metric as 

4(A) EE {d,^d,<t>), ~ {d.<p)s {d,<P)s . (108) 



IV. APPLICATION TO SIGNALS FROM 
ISOLATED NEUTRON STARS 

In the following we apply the general framework of the 
previous sections to a particular class of gravitational- 
wave signals, namely GWs emitted from isolated spin- 
ning neutron stars. Restricting ourselves to a specific 
signal model allows us to explicitly compute the metric 
and evaluate different approximations for various cases 
of interest. It also allows us to compare the metric pre- 
dictions to measured mismatches using a realistic search 
code for computing the J-"-statistic. 



A. The phase model 

GWs emitted from isolated spinning neutron stars can 
be described by a very simple phase model, namely 



k=0 ^ ^ >■ 



(109) 
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where /^'^•' is the fc-th time derivative of the intrinsic sig- 
nal frequency /(t) in the solar-system barycenter (SSB), 
and T-^(i) is the arrival time in the SSB (with respect 
to the reference time Trof) of a wavefront reaching the 
detector X at time t. Neglecting relativistic effects, this 
relation is simply given by 



(110) 



where r^{t) is the position of detector X with respect to 
the SSB. This can be separated into an orbital and a spin 
component, namely 



rh(t) 



(111) 



where rovh is the position of the earth in the SSB, and 
rf^ijj is the position of the detector on earth. The unit 
vector n denotes the sky position of the source, which 
can be written as 



(cos (5 cos a, cos (5 sin a, sin ^) 



(112) 



in terms of the equatorial coordinates right ascension a 
and declination S. For this phase model we have the 
derivatives 



90^(t) „ 



(113) 



dn 



27r 



{k + iy. ' 



k=0 



The intrinsic signal frequency fir) can usually be as- 
sumed to be a slowly varying function of time, and typ- 
ically one only needs to include a small number of spin- 
down coefficients /^'^^ between s — and s — 3, say. 
With these phase derivatives and the explicit antenna- 
pattern functions a{t), b{t) (Eqs.(12),(13) in 0), we 
can numerically compute the different metrics derived in 

sect.irm 

The corresponding time integrals involved in these 
expressions were computed numerically using a Gauss- 
Legendre quadrature of order 2000. This results in a (sat- 
urated) convergence of the metric components gtj and the 
antenna-pattern integrals A, B,C a.t the level of a relative 
precision of 10"^*^, while the corresponding mismatches, 
m — AA*AA^, only converge to a level of about 10~^. 
The weaker convergence of m is can be attributed to the 
fact that the metrics are highly ill-conditioned matrices, 
as will be discussed further in Sect. IIV E 21 



B. Natural units for the Doppler offsets AA 

We can easily find the scaling and the order of magni- 
tudes of the metric components: keeping only the dom- 
inant scaling-terms in the phase derivatives (|113l) . 1114(1 
we find 



df{k)(t)- 



dn4'^ frorh{t)/c. 



(115) 



Taylor-expanding rorb(i) for "short" observation times 
T <C 1 year, and neglecting the antenna-pattern func- 
tions a{t) ~ h{t) ^ 0(1), we find the following scaling for 
the metric components (I107|l : 



{fV/crr^e^e,, 
T'^+i fTV/c e, , 



(116) 



where e' is the unit vector along the orbital velocity V^. 
These estimates allow us to read off the (minimal) scaling 
of the metric determinant g = det gij and of the corre- 
sponding volume measure y^, namely 



fc+i 



(117) 



fe=0 



where s is the number of spindown coefficients to in- 
clude. Note that the dominant term Hllt)|) of the sky 
metric 5„i„j is degenerate and will only be regularized 
by higher order terms. Therefore the scaling of the num- 
ber of sky templates will be at least of the order 0(T^), 
but is likely to be higher. If we consider the case of a 
single spindown, i.e. s — 1, we find the scaling for the 
volume (and therefore number of templates) of the order 

In order for these metric components to be dimension- 
less and of order unity, we should measure the Doppler 
offsets AA^ in terms of "natural units" AAg, namely 



A/, 



(k) 



l/T'+\ Ano = c/ifTV), 



(118) 



where V/c « 10 ^. We denote the Doppler offsets in 
natural units as 



AA' = AAVAAj, . 



(119) 



Note that n is the unit vector pointing to a sky location, 
and using its expression ((112|l in terms of equatorial co- 
ordinates a, 6, we find 



|An|2 = cos^^ (Aa)2 + {ASf 



(120) 



For small offsets, |An| is simply the angular distance cor- 
responding to offsets Aa and Ai5. We therefore define the 
natural-unit angular offset Afi on the sky as 



An = 



\An\ 
Ano 



(121) 



In these units, the statement of "small" Doppler offsets in 
the expansion (|43l) is really meaningful, i.e. the validity 
of the metric approximation is guaranteed in the regime 



A/C') < 1 , Ah-^1 



(122) 



However, as will be seen later, this condition is in many 
cases not necessary for the approximate validity of the 
metric. Note that the metric iso-mismatch ellipses on 
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the sky are highly anisotropic (the dominant contribu- 
tion shown in H116|) is in fact degenerate), and the actual 
metric scale can therefore deviate significantly from Ail, 
depending on the direction of the angular offset. In fact, 
An corresponds to the larger of the two angular eigenval- 
ues on the sky (resulting in the smaller mismatch-scale), 
but it nevertheless captures the correct scaling of the sky 
metric. 



C. Numerical exploration of the parameter space 

The results in the following sections are based on 
Monte-Carlo simulations of the parameter space of sig- 
nals and Doppler offsets. Here we summarize how the un- 
derlying random distributions are generated. For the sig- 
nal amplitudes A'^ , given in , we pick cos l and ip with 
uniform distributions from the ranges cos l € [~ 1 , 1] and 
ip e [0, tt]. We simply fix the overall amplitude to /iq = 1 
and the initial phase to 0o — 0, as these parameters 
have no effect on the metric. The Doppler parameters 
A* of the signal are picked (with uniform distributions) 
from the ranges / e [100,200] Hz, / e [-10-^10-9] s~^, 
a € [0,27r] and sin (5 G [—1, 1]. Note that S is picked in 
such a way as to obtain an isotropic distribution of points 
{a, (5} on the sky sphere. 

The appropriate selection of random Doppler offsets 
AA* is more delicate, because their distribution should 
be roughly "isotropic" in some metric sense, and the cor- 
responding J^-mismatches should lie in a "reasonable" 
range, where the metric approximation is applicable, e.g. 
mj^ < 0.5 say. Note that the second requirement is ir- 
relevant for the purpose of studying intrinsic properties 
of the metric and for comparing different metric approx- 
imations with each other. However, it is essential for 
a comparison of the metric predictions to measured J-- 
mismatches, which will be described in Sect. IIV Ul 

We use the following algorithm to generate suitable 
Doppler offsets AA*, satisfying the above requirements: 

1. Pick the signal parameters and A* as described 
above, and compute the corresponding J^-metric 
gfj{cos L, Ip; A) at this parameter-space point. 

2. Pick a random offset vector A A * in natural units 
with a uniform distribution in the hypercube, i.e. 
AA * G [—1,1]. The corresponding dimensional 
offsets A A'* are therefore uniformly distributed in 
A/' e [-A/o,A/o], {Aa',A<5'} e [-Ano,Ano], 
and A/'e [-A/o,A/o] 



Normalize AA'* using the metric, i.e. construct 
e\ = AA'VIAA'I, where |AA'| = 



.g^AA'^AA'^' 



4. Pick a distance d = ^Jm with uniform distribution 
in d e [0.1, Vo^]- 

5. The resulting Doppler offset is then obtained as 
AA* = d e\. 



The lower bound of mj: > 0.01 on the generated JF- 
mismatches is chosen in order to avoid numerical errors 
dominating the relative differences between mismatches, 
which happens especially when comparing to measured 
mismatches. We typically generate about 250,000 ran- 
dom choices of signal parameters and offsets, for obser- 
vation times ranging from T = 12 hours to T = 200 hours 
in steps of 4 hours, which corresponds to about 5000 tri- 
als for each choice of observation time. By construction, 
the distribution oi mjr is such that ^mjr is uniformly 
distributed in the range [0.1, -\/0.5]. The distribution of 
the Doppler offsets AA in natural units is found to be 
approximately Gaussian with zero mean and a standard 
deviation of about 0.3. There are, however, a few percent 
of offsets reaching up to a few hundreds in natural units, 
which corresponds to highly degenerate directions of the 
metric. 



D. Intrinsic uncertainty of the JT- metric 

As discussed in Sect. Ill Al the reduced parameter space 
of the J^-statistic is the Doppler space A, while the am- 
plitude space A has been "projected out" by maximiza- 
tion. For the problem of covering the Doppler-parameter 
space with templates, the amplitudes represent unknown 
external parameters, and so we cannot directly use the 
J^-metric family gf^{A, A) for the covering problem. How- 



ever, as shown in Sect. 



we can compute an average 



J-"-metric 5^(A), and a corresponding intrinsic relative 
uncertainty CT:r(A, AA), which are both independent of 
the unknown amplitudes A. This uncertainty range in 
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FIG. 1: Iso-mismatch (m = 0.1) ellipses of the jF-metric fam- 
ily gfj{cos X) for 100 random choices {cost, and of 
the average metric 'gfj{X). Parameter-space cut along the sky 
plane {a, 5}. The top and right-hand axes show the corre- 
sponding offsets in natural units. [Parameters: / = 100 Hz, 
a = 1.45 rad, 5 = rad, / = 0, detector = 'LI', GPS start- 
time to = 810720013 s, duration T = 50 hours.] 
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-3 3 
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FIG. 2: Same as Fig. for a parameter-space cut along the 
{/, /} plane. The top and right-hand axes show the corre- 
sponding offsets in natural units. 



the predicted mismatch is illustrated in Fig. ^ and Fig.|21 
which show different parameter-space cuts through the 
iso-mismatch hyper-ellipsoids for the average metric gfj 
and for 100 randomly picked members of the J^-metric 
family (/-^(cos l, Note that the intersection of all these 
mismatch hyper-ellipsoids myr < 0.1 corresponds to the 
worst-case mismatch region m^'^^ < 0.1, which is clearly 
not a hyper-ellipsoid. Fig.Elshows the distribution of the 
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FIG. 3: Distribution of the relative intrinsic uncertainty 
of the jT-mismatch as a function of observation time T, for 
the single-detector case. Plotted are the mean, standard de- 
viation and extremal values of the distribution for each value 
of T. 

relative uncertainty CTjf(A, AA), as a function of the ob- 
servation time T for randomly chosen signal parameters 
and offsets (using the algorithm described in Sect. lIVC)l . 
These results are for the single-detector case only, the 
corresponding dependence on the number of detectors is 
investigated in Sect. lIVFl We see in Fig.|21that for short 



observation times, of the order of T 12 hours, the 
relative uncertainty can be as large as 100%. For longer 
observation times of the order of a few days, crjr decreases 
substantially, but even for T ~ 8 days, the intrinsic un- 
certainty can still reach up to about 20 %. 



E. Comparing different metric approximations 

For the following comparisons of different metric ap- 
proximations, it will be useful to define the relative error 
eia, b) between two quantities a and b as 



e(a, b) 



(|a| + |6|)/2 



(123) 



which is the harmonic mean of {a — b)/\a\ and {a—b)/\b\. 
This definition has the advantage of being bound within 
[—2,2] even for large deviations a ^ b or b ^ a, while 
it agrees with the more common definitions of relative 
errors for small values |e| ^ 1. 



1. Monte-Carlo study of mismatches 

As we have seen in Fig. |31 the intrinsic relative uncer- 
tainty (TjF decreases substantially on the timescale of a 
few days, corresponding to a convergence gfy The 

average metric 'gf^ is therefore expected to be an increas- 
ingly reliable approximation of the full J-"-metric with 
longer observation times T. This is indeed the case, as 
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FIG. 4: Distribution of relative errors e{mr,rny^) of the 
(single-detector) average .7^-mismatch mjr with respect to the 
exact .7-"-mismatch mjr as a function of observation time T. 
Plotted are the mean, standard deviation and extremal values 
of the distribution for each value of T. 

can be seen in Fig. 0] which shows the distribution of rel- 
ative errors e(?7ijr, mjr) as a function of observation time 
T. This distribution was obtained using random choices 
of the signal parameters and Doppler offsets as described 
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in Sect. IIVO The standard deviation of the relative er- 
rors is about 20 % for short observation times T <1 day, 
but rapidly decreases to about 1 % for timescales of a few 
days. Note that these errors are somewhat smaller than 
could have been expected from the intrinsic uncertainty 
CTjF shown in Fig.|3 which indicates that mismatches near 



the extremal values m 
closer to the average rrijr. 



mm max 



are less likely than those 



As shown in Sect. 1111 El the JF- metric family con- 
verges to the orbital metric gf^ for very long observa- 
tion times r 3> 1 day, but it is not obvious on which 
timescales this convergence happens in practice. Fig. |31 
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FIG. 5: Distribution of relative errors e(m;F,morb) of the or- 
bital metric with respect to the (single-detector) .7^-metric 
as a function of observation time T . Plotted are the mean, 
standard deviation and extremal values of the distribution for 
each value of T. 



shows the relative errors eijrijr, morb) of the orbital met- 
ric compared to the exact J^-metric, as a function of ob- 
servation time. Note that a plot with the distribution of 
e{mjr^m^) would look virtually indistinguishable. These 
two metric approximations perform nearly identically as 
far as predicting mismatches is concerned, as can be 
seen in Fig. showing the distribution of relative errors 
e(m0,morb)- Note that the orbital as well as the phase 
metric show a tendency to overestimate the mismatch for 
short observation times T < I days, which is apparent 
in Fig.|Sl Summarizing, we conclude that the "mismatch 
quality" of the orbital and of the phase metric seem vir- 
tually identical, and both approximations perform only 
sightly worse than the average J^-metric shown in Fig.^] 
However, as will become clear in the next section, there 
is an important difference between the phase metric and 
the orbital metric, which does not manifest itself in this 
Monte-Carlo study due to the smallness of the relevant 
parameter space in which the difference becomes appar- 
ent. 
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FIG. 6: Distribution of relative errors e{m^,morb) of the or- 
bital metric with respect to the phase metric as a function of 
observation time T. Plotted are the mean, standard deviation 
and extremal values of the distribution for each value of T. 



2. Metric determinants and eigenvalues 

A complementary way of comparing different metrics is 
to look at their eigenvalues and determinants. These in- 
variant quantities characterize the semi-major axes and 
the volume of iso- mismatch hyperellipsoids, which are 
important properties for the template covering of the pa- 
rameter space. We denote the four metric eigenvalues as 
5i ^ 92 > 33 > .94- The metric for isolated neutron- 
star signals is quite generally represented by highly ill- 
conditioned matrices in the default parameter-space co- 
ordinates (a, (5, /, /, ...). That means that gi is typically 
many orders of magnitude larger than 54, correspond- 
ing to very thin and elongated mismatch hyperellipsoids. 
This property can be characterized by the so-called condi- 
tion number, defined as k = gi/g^. Well-conditioned ma- 
trices have K ~ 0(1), while the metrics encountered here 
typically have K.{gfj) ^ 10^° in SI units, and i^igf^) ^ 10^ 
in natural units (cf. 



Sect. I1VB(I . Ill-conditioned matri- 
ces can strongly amplify numerical errors in computing 
their maps, inverses, determinants or eigenvalues. One 
has to be very careful when handling such matrices nu- 
merically, as the results can be quite unreliable. Using 
SI units, it turns out to be all but impossible to com- 
pute the determinant or the eigenvalues using standard 
double-precision arithmetic, which suffers from complete 
loss of significant digits in this case. In natural units, 
however, these quantities can be computed (albeit not 
with good precision, and not for all points in parame- 
ter space), and the results for one example are shown in 
Fig. 13 and Fig. |S1 Contrary to the apparent fast con- 
vergence of the orbital- and the ^-mismatches shown in 
Fig. [SJ we see a very different picture for the metric de- 
terminants in Fig. \7\ For timescales of a few days, the 
determinant of the orbital metric is a few orders of mag- 
nitude smaller than the determinants of the other met- 
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FIG. 7: Determinants of different metric approximations (in 
natural units) as functions of observation time T. [Parame- 
ters: / = 100 Hz, / = 0, a = 1.0 rad, 5 = 0.5 rad, detector = 
'HI', GPS start-time to = 792576013 s] 



ric approximations, and the expected convergence H1U()|) 
takes place only for timescales longer than T > 1 month. 

How can this be reconciled with the apparently much 
faster convergence of the mismatches in Fig. (Sp* In order 
to better understand this, let us look at the four eigen- 
values gi as functions of T, which is shown in Fig.|Sl We 
see that the largest discrepancy of the orbital metric oc- 
curs for the smallest eigenvalue, (74, corresponding to the 
most degenerate direction of the metric. The largest two 
eigenvalues gi , 32 agree well, and 33 only differs by a fac- 
tor of a few. In order for the mismatch to be affected by 
the smallest eigenvalue, we would have to pick a Doppler 
offset that is very closely aligned with the most degen- 
erate principal axis, as any appreciable offset along the 
other axes would easily dominate the total mismatch. A 
very rough order-of-magnitude estimate of the probabil- 
ity p of picking such a direction yields an upper bound 
of p < \/(54/5i)(54/52)(ff4/ff3) < 10"*^. It is therefore 
very unlikely to pick a Doppler offset for which the mis- 
match is dominated by the smallest eigenvalue. This is 
consistent with the fact that in about ~ 10^ trials we 
did not see any cases in which the orbital metric had 
dramatically underestimated the mismatch, i.e. where 
e(m;(r,TOorb) 2 in Fig. 

It is interesting to note that while the phase metric 
seems virtually identical to the orbital metric for "al- 
most all" directions in parameter space, its determinant 
(and eigenvalues) agree much better with the (average) 
^-metric. The small effect of the spin motion of the earth 
is negligible for most directions in parameter space, ex- 
cept for the most degenerate one, where it substantially 
reduces the degeneracy. 



F. The multi-detector JT- metric 

As discussed in Sect. IIII CI the parameter-space reso- 
lution of the multi-detector ^-metric does not scale with 
the number of detectors. Instead, the effect of combining 
detectors coherently results in a noise-weighted average 
of contributions from different detectors. This averaging 
operation H59() would be expected to decrease the effects 
of the antenna-pattern functions a?^{t) and b^{t), as well 
as the detector-specific Doppler modulation A(fP^{t) of 
the signal phase H96|) . In Fig.|51we see indeed that both 
the intrinsic uncertainty aj^ of the J^-metric family, as 
well as its relative difference to the orbital metric de- 
crease with the number of detectors. These results are 
based on a Monte-Carlo simulation with ^ 40, 000 ran- 
domly chosen parameters (see Sect. IIV Cll . for a fixed 
observation time of T = 55 hours, and using between 
one and six coherently combined detectors. For the sake 
of this example, we have made the (obviously unrealis- 
tic) assumption that all 6 detectors have the same noise 
floor, otherwise the convergence would be much weaker. 
Note that we would find exactly the same mismatch con- 
vergence for the phase metric as that shown in Fig. El^b) 
for the orbital metric. 



G. Comparing metric predictions to measured 
mismatches 

In order to validate the theoretical JF-statistic mis- 
match mjr, derived in Sect. IIII Bl we compare it to 
the measured relative SNR loss mo in a simulated mis- 
matched search. This is done by first generating (us- 
ing lalapps_Makef akedata 15j) a fake signal with pa- 
rameters {A, As} picked at random (see Sect. IIV^ . We 
then measure (using lalapps_ComputeFStatistic [isj ]') 
the perfectly-matched SNR p(0) at the signal location 
As, and the mismatched SNR /o(AA) at an offset Doppler 
position As -I- A A. The offsets A A were picked at ran- 
dom using the algorithm described in Sect. IIV CI The 
measured J^-mismatch mo is then given simply by the 
definition Eq. (|46|l . The comparison of the measured 
mismatches mo to the theoretical approximation reveals 
a problem with the metric approximation for large angu- 
lar offsets (in natural units), in particular for Af2 > 5. 
The origin of this "metric failure" can be understood in 
terms of the metric curvature and will be discussed in 
more detail in IIV C 21 In the meantime we simply re- 
move this known source of errors by excluding points 
with Ail > 5, which affects less than 1% of the Monte- 
Carlo trials. Fig. 1101 shows the measured mismatches mo 
as a function of the prediction mjr, for observation times 
of (a) T = 12 hours and (b) T = 60 hours. We see that 
there is a substantial scattering of measured mismatches 
in the case of T = 12 hours, which has virtually disap- 
peared for r = 60 hours. This second type of "metric 
failure" , which only affects short observation times, will 
be discussed in IIV G ll 
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FIG. 8: Square roots of the four eigenvalues gi of different metric approximations (in natural units) as functions of observation 
time T. [Same parameters as in Fig. |7| 
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FIG. 9: Distribution of (a) intrinsic relative .7^-metric uncertainty ay^ and (b) relative errors e{mr,mc,ih) as functions of the 
coherent combination of (assumed equal-noise) detectors, for T = 55 hours. Plotted are the mean, standard deviation and 
extremal values of the distribution for each detector combination. The assumed detectors (location and orientation only) are: 
'HI' = LIGO Hanford, 'LI' = LIGO Livingston, 'Gl' = GEO600, 'VI' = Virgo, 'Tl' = TAMA, 'PI' = Cahech 40m. 



Another effect seen in these figures is a systematic de- 
viation of rriQ with respect to with increasing mis- 
matches, which becomes noticeable at around my^ > 
0.15. This deviation would be suspected due to higher- 



order corrections 0(AA^) with respect to the local met- 
ric expansion, and is found to be roughly independent 
of the observation time T. We can approximate this sys- 
tematic deviation by an empirical quadratic correction of 
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FIG. 10: Measured mismatches mo versus predicted .7-"-mismatches niT for observations times of 
T = 60 hours, omitting points with Afi > 5. The dashed Une indicates the identity function. 
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the form 



mjr^qn{mjr) = myr - 0.38 m 



(124) 



which can be used in comparisons to measured mis- 
matches, as it compensates for the systematic "drift" in 
mp. Note that this deviation makes the metric a conser- 
vative over-estimate, as the actual mismatches tend to 
be smaller than the predicted ones. 



1. Metric problems for T < 1 day 

The strong scattering of mismatches seen in Fig- lior a) 
for short observation times can be attributed to an in- 
trinsic property of the J^-statistic: namely, the "local" 
parameter-space structure for short observation times 
is not very well approximated by a quadratic form in 
offsets. Note that obviously the metric approximation 
is valid in the strict local sense of sufficiently small 
offsets, but in practice we are more interested in a 
finite "local" region of small mismatches, m < 0.5 
say. An illustrative example for this is shown in 
Figs, [m lEl where we have measured the J^-statistic 
(using lalapps_ComputeFStatistic) over a ~ 10~^ Hz 
band around the true signal frequency, with the other 
Doppler parameters {a,S,f) held fixed at their correct 
values. For an observation time of T = 5 hours, the 
quadratic decrease predicted by the JF-metric is clearly 
not a good approximation, as can be seen in Fig. 1111 
depending on the sign of the frequency offset, the met- 
ric would either substantially under- or overestimate the 
mismatch. This effect decreases rapidly with observation 
time, and for T = 30 hours (and the same signal), the de- 
crease of the J^-statistic is "locally" much better approx- 
imated by the metric, as seen in Fig. E| This effect can 
be understood as follows: the quadratic decrease of the 
J^-statistic predicted by the metric has a typical width in 




FIG. 11: Comparison of the measured decrease in .7^-statistic 
(CPS) as a function of frequency offset A/, and the corre- 
sponding metric prediction 1 — mjr, for an observation time 
of T = 5 hours. [Parameters: / — 100 Hz, cos t = 0.95, 
ijj = 2.25, a = 5.98, S = 0.09, / = 0] 



frequency of A/o ^ 1/T (cf. Sect. IIV S^ . However, the 
spin motion of the earth creates "side-lobes" in frequency 
space at offsets of a few /spin = l/day~ 10~^ Hz, which 
are clearly visible in Fig. 1121 For observation times much 
smaller than a day, these side-lobes are therefore not well 
separated from the main peak and will substantially al- 
ter its quadratic form, e.g. by creating "plateaus" and 
steeper "cliffs" , as seen in Fig. ^2 As the observation 
time gets longer, i.e. T > 1 day, the main peak becomes 
well separated, and therefore well described by the metric 
approximation. 
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FIG. 12: Same as Fig. 1111 for an observation time of T 
30 hours. 



2. Metric problems for AQ > 5 

As already mentioned, another problem with the met- 
ric approximation affects points with large angular off- 
sets Ai7, and is present even for long observation times. 
This can be seen in Fig. 1131 showing the relative errors 
e{mo I w.F,qu) as a function of angular offset Ail, for ob- 
servation times T > 48 hours. For angular offsets Ail < 
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FIG. 13: Relative errors e(mo,mjr qu) versus angular offset 
AQ in natural units, for observation times T > 48 hours. For 
better readability of the figure, the displayed range of angular 
offsets is limited to Af2 < 40. 



tion. This curvature is closely related to the global "cir- 
cles in the sky" (CiS) structure discussed in Q, as the 
metric ellipses on the sky are "tangential" to these cir- 
cles. To first order, the CiS are described by the equation 



/I 



V/c)=fJl 



V/c 



(125) 



where / and n are the "target" frequency and sky posi- 
tion respectively, while fs and rig are the corresponding 
signal parameters. The problem stems from the curva- 
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FIG. 14: Metric "failure" for large angular ofi'set AO, = 84 
at r = 68 hours. The point 'x' indicates the signal location, 
shows an offset location on the m = 0.05 iso-mismatch 
ellipse with a measured mismatch of mo ~ 1. The dashed line 
indicates the "circle in the sky" l|125|l defined by the signal 
location. The top and right-hand axes show the offsets in 
natural units. [Parameters: / — 186.34 Hz, a — 3.93, 5 — 
-0.33, / = 1.7 X 10-^°s-^ cost = 0.4, tp = 1.65 ] 

ture of the JF-statistic circles in the sky-coordinates a, 6, 
which is why the CiS is only locally well approximated 
by the respective mismatch ellipses. This effect is seen 
clearly in Fig. 1141 which shows an extreme example of 
such a metric "failure" due to large Afi. In this fig- 
ure, indicates a sky position on the mjr = 0.05 iso- 
mismatch ellipse, which has an angular offset from the 
signal of Afi — 84. Contrary to the predicted mismatch, 
the measured mismatch at this point is mo ~ 1 , as the !F- 
statistic decreases rapidly away from the CiS (indicated 
by the dashed line). 



5, the relative errors stciy below si^uiQ^TTij^' qn ) < 0.1, but 
with increasing Afl, the errors start to spread out sub- 
stantially. 

The reason for this metric failure can be traced to the 
curvature of the metric on the sky. While the metric el- 
lipses have constant orientation as functions of frequency 
and spindown, their orientation changes with sky posi- 



3. Summary of metric validation 

Fig. II 5l summarizes the relative errors s{mo, mj^.qu) be- 
tween the (quadratically corrected) predictions (|124|l and 
the measured mismatches mg, as a function of observa- 
tion time T. We see that, as expected, excluding large 
angular offsets (Afi > 5) yields a substantially improved 
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FIG. 15: Distribution of relative errors e(mo 
function of observation time T. Plotted are the mean, stan- 
dard deviation and extremal values of the distribution for 
each value of T, excluding points with AQ > 5. The points 
'x' and '+' indicate the respective maximum and minimum 
errors when including results with An > 5. 



agreement between the metric prediction and the mea- 
surements, as it ehminates the type of "metric failures" 
discussed in IIV G 21 We also see that the relative er- 
rors can still be quite large for short observation times 
T < 1 day, as discussed in IIV G II and that these er- 
rors decrease rapidly with T. For observation times of 
T > 1 day, the average error is below a few percent. 
We can therefore conclude that the agreement of the J-- 
metric with the measurements is very good in the domain 
of applicability of the metric approximation. 



V. DISCUSSION 

We have derived a formalism for the general 
parameter-space metric of the multi-detector J^-statistic, 
and we have explicitly computed the metric for signals 
from isolated spinning neutron stars. We have shown 
that there exists a family of ^-metrics, parametrized by 
the two (unknown) amplitude parameters V' and cos i. 
We explicitly derived the extremal "mismatch bounds" 
(i.e. the maximum and minimum possible mismatches) 
of the J-"-metric family, and wc introduced an average JF- 
metric, which is independent of the unknown amplitude 
parameters. We have shown that the multi-detector J-- 
metric does not scale with the number of detectors. Com- 
bining detectors coherently therefore does not increase 
the required number of templates. In the long-duration 
limit (T > 1 month), wc found that the JF-metric fam- 
ily converges towards a simple orbital metric 5°^*^, which 
neglects both the amplitude modulation and the phase 
modulation caused by the diurnal rotation of the Earth. 

Both the orbital and the closely related phase met- 
ric provide relatively good mismatch approximations in 



practice, and the quality of these approximations im- 
proves with longer observation times and with the num- 
ber of coherently-combined detectors. 

The orbital metric, however, while virtually identical 
to the phase metric for almost all directions in parameter 
space, was found to be substantially more degenerate for 
observation times shorter than a month. This has impor- 
tant consequences for the covering problem of the param- 
eter space and requires further study. Finally, we have 
identified two regimes in which the local metric approxi- 
mation itself is not very reliable: namely, for observation 
times T < 1 day, and for large angular offsets Afl > 5 
(in natural units). 



APPENDIX A: ALTERNATIVE DERIVATION OF 
THE JP-METRIC 

A more elegant derivation of the ^-metric (|47|l can 
be obtained j20| | by projecting the full parameter-space 
metric in {A, A} into the reduced parameter space A of 
the JF-statistic. We denote the full parameter space as 
9 — {A, A}, and we use index conventions 9°- — {A'^, A*}. 
Recall the form of the log-likelihood ratio IjlSfl . namely 

\nA{x;9)^{x\s{9))-l{s{e)\s{9)) . (Al) 

If the data contains a signal with parameters Og, i.e. 
x{t) = n{t) + s{t;9s), and if the target position 9 is 
"close" to the signal location, i.e. 9 = 9s + A9 for "small" 
A9, then the expectation value of In A can be expanded 



2E[lnAi9\9s)] = is\s)~idas\dbs) A9'' A9'' + 0(3) , (A2) 
so the full parameter-space metric gat is found as 



gab 



{daS\dbS) 

{s\s) 



(A3) 



This expression is sometimes referred to as the normal- 
ized Fisher matrix . As mentioned in the introduction, 
this "canonical" metric differs from a definition often 
found in the literature (e.g. T, 2,'T2]), which is based on a 
somewhat more ad-hoc measure of the "match" , namely 
M = {s{9,)\s{9)), instead of the full log-likelihood l|IT|l . 
As it turns out, both definitions result in the same phase 
metric (|107f) when considering constant-amplitude sig- 
nals (after minimizing the mismatch (|A3(I over the un- 
known amplitude). In general, however, the canonical 
definition ljA3|l . and correspondingly (|46() . is more di- 
rectly relevant to the covering problem, as it describes 
the relative loss of detection statistic. For the assumed 
general form of the signal we have 



is\s) = A^M^^.A" 



(A4) 
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and the respective derivatives with respect to the ampli- 
tude and Doppler subspaces are given by 



ds 



(A5) 
(A6) 



The full metric gab therefore consists of the three blocks 
(with respect to the two subspaces): 



9at - p-'{0) 



(A7) 



in terms of hij and TZi defined in (|40|l and (|31|l . The 
^-metric gf^ in the Doppler subspace can be regarded 
as the distance corresponding to given Doppler offsets 
A9^ — AA*, minimized over the amplitude offsets AQf^ = 
AA^', i.e. 



gf^AV AX^ = min^,. gab AO" Ae" 



(A8) 



This can be minimized trivially, since it is a quadratic 
function in AA'^, and the "compensating" amplitude 
mismatches are 



Ae'^ = -r''gu:A0\ 



(A9) 



where g'^''' is the inverse matrix of g^^, i.e. g'^°'gaiy — S^- 
Inserting this into ljA8|l we obtain 



(AlO) 



which corresponds to the projection of the full metric gab 
into the Doppler subspace. Using the explicit compo- 
nents (jA7|) . we find 

gf^iA) = p-^iO)A[h,, ~n^,M-^n,]A (All) 



A-M-A 



(A12) 



in perfect agreement with the earlier result (|48|l . which 
was obtained in a more straightforward, but somewhat 
more tedious calculation. 



APPENDIX B: .7^-METRIC FOR 
LOW-FREQUENCY SIGNALS 

In the case of low frequencies /g (which would be rele- 
vant for LISA) and/or short observation times T, where 
fsT 10^, we cannot use the simplifying approxima- 
tion of Sect. IIII Dl We can nevertheless proceed in the 
same way: using the expansion 1)64(1 and keeping only 
leading-order terms, we can express (|40ll as 
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in terms of the four independent components 

{diadja)g + {a? di(l)dj(f>)g , 
(dibdjb)g + (6^ di(t)dj<t))g , 



P.2 



i^^- = {diadjb)g + {abdi(l)dj 



(B2) 



fj = {bdiadj(l))g - {adibdjc, 



with implicit symmetrization in i, j. In a similar manner 
we calculate TZ^i, i , defined in ((^{1(1 , which yields 
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S-^T 
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where the 2x2 matrices TZi and TZi are defined as 
TZ^ E 



pll J^12 

■21 7j22 ) 5 and 7?.i 



7j2i 7J2 



PP R}^ 

nl4 n24 
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in terms of 





(a 5ia)g 


pl2 


= {adib)s , 


d21 

— 




p22 






[a'^di(j))g 


Ri 


= (629,0)s 






{a b di4> 
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(B5) 



Note that Ri only contains derivatives of the antenna- 
pattern functions a^, b^ (which where neglected in 
Sect. IIII Dll . while Ri only contains derivatives of the 
phase. As a consequence of this block structure, one 
finds 



i'^TZ^M-^TZA ^-S-^T 



^ Qij Qij 
Qij Qij 

4. 











with the four independent components 
DQl^ = A [PfPf + ^14^14] ^ 5 [^11^11 
-2C[R}^Rf+RfR]^] , 



i?Q2^ = A[i?f4^24 



i?f2i?22] ^ 5 [^14^14 





% ' 
(B6) 

- Rl'R]'] 



R^'R]'] 



d12 Dill 



-2C [Rl^Rf + PpPf ] , 
- A [Rl^Rf + i?f i?f ] + B [Rl^Rf 

n r 7?24 E?13 I d14 nl4 i nil n22 i pl2 7j211 

— O [rtj rCj -h Kj -|- Kj ii^- -|- rCj rC^ J , 

Q4. ^ ^ [^21^24 _ ^14^22] + 5 [^11^14 „ ^12^13] 

r'r??ll??24 Rl3 p22 I pl4 p21 pl2 pl41 

— O [rCj rCj — Kj rt^- -h tx^ Hj — rCj ri^ J . 

(B7) 

Putting all the pieces together, we find Qij , defined in 
(05, in the form: 



Guv ~ 7:5 
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— m 
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—m 




m 


I 


4 

rn 





m 


m 



(Bl) in terms of the /owr components (r = 1, 2, 3, 4): 
m'- ^ (Pr - Q-^ AA'AA^ . 



(B8) 



(B9) 
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As discussed in Sect. IIII( the extrema fhjr are the eigen- 
values of Ai^^ ■ Q, which are the solutions of 

= detTW"^ dct[g -fhjrM] , (BIO) 

or equivalently 

— {m^ — rhjrA){rn'^ — fhjrB) 

-(m3-m^C)2-(m4)2. (Bll) 

We see that there are again maximally two independent 
eigenvalues, namely 



which is formally very similar to the earlier result (|90|l . 
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